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Using large-scale Monte Carlo computations, we study two versions of a (1 + 1)D Z4-symmetric 
model with ohmic bond dissipation. In one of these versions, the variables are restricted to the 
interval [0,27r), while the domain is unrestricted in the other version. The compact model features 
a completely ordered phase with a broken Z4 symmetry and a disordered phase, separated by a 
critical line. The noncompact model features three phases. In addition to the two phases exhibited 
by the compact model, there is also an intermediate phase with isotropic quasi-long-range order. 
We calculate the dynamical critical exponent z along the critical lines of both models to see if the 
compactness of the variable is relevant to the critical scaling between space and imaginary time. 
There appears to be no difference between the two models in that respect, and we find 2 « 1 for 
the single phase transition in the compact model as well as for both transitions in the noncompact 
model. 
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I. INTRODUCTION 

The standard way of introducing dissipation in a quan- 
tum mechanical system is to couple some variable de- 
scribing the system to the degrees of freedom of an ex- 
ternal environment^. The environment is modeled as a 
bath of harmonic oscillators which couple linearly to the 
system variables. The oscillator degrees of freedom, ap- 
pearing in the action to second order, may be integrated 
out to produce an effective theory for the composite sys- 
tem given in terms of the system variables. 

The presence of a dissipative term introduces strongly 
retarded (nonlocal in time) self-interactions of the sys- 
tem variables. This long-range interaction in imaginary 
time may have serious consequences for the quantum crit- 
ical behavior of the system. This effect can usually be 
described by a dynamical critical exponent z defined by 
the anisotropy of the divergence of the correlation lengths 
at criticality, ~ C^i where ^ and are the correlation 
lengths in space and imaginary time, respectively. An 
Ising spin chain with site dissipation was shown by ex- 
tensive Monte Carlo simulations in Ref. 2, to have z ~2. 
The same model, augmented to two spatial dimensions, 
was investigated by the present authors in Ref. [31 The 
result z = 1.97(3) suggests that the dynamical critical 
exponent is independent of the number of spatial dimen- 
sions, in agreement with naive scaling arguments which 
make no reference to dimensionality.!^ On the other hand, 
when coupling the reservoir to bond variables involving 
Ising spins the dissipation term was found to be irrele- 
vant to the universality class, i.e., z « 1."^^ In general, 
dissipation suppresses certain types of quantum fluctu- 
ations, though the larger value of z for site dissipation 
signifies that bond dissipation is far less effective than 
site dissipation in reducing fluctuations. 

Ohmic dissipation in terms of gradients or bonds is 
common in models describing shunted Josephson junc- 
tions or granular superconductor systems .^SEj Here, the 



bonds represent the difference of the quantum phases 
between the superconducting grains. In this context it 
is well-known that the coupling of the environment to 
the system may affect the natural domain of the system 
variables For Josephson junctions, this means that the 
domain of the phase variables reflects quantization of the 
charges on each superconducting grain. If the charges are 
quantized in units of Cooper pairs, 2e, the domain of the 
quantum phase is 27r-periodic. Ohmic shunting leads to 
an unbounded —00 < 6 < 00 domain,-^ reflecting a con- 
tinuous transfer of charges across the junction. We will 
henceforth refer to the variable defined on a restricted 
27r interval as compact and the extended variable as non- 
compact. 

Moreover, dissipation in terms of bonds has also been 
proposed in an effective model describing the low-energy 
physics of fiuctuating loop currents to describe anoma- 
lous normal state properties of high- Tc cuprates.^^ A 
quantum statistical mechanical model for such degrees 
of freedom has been derived from a microscopic three- 
band model of the cuprates.l^ The classical part of the 
derived action consists in its original form of two species 
of Ising variables within each unit cell, coupled by a four- 
spin Ashkin- Teller term. This model has been proven, 
through large-scale Monte Carlo simulations, to support 
a phase transition with a non-divergent non-analyticity 
in the specific heat on top of an innocuous background. 
The breaking of the Ising-like symmetry describes a sug- 
gested ordering of loop currents upon entering the pseu- 
dogap phase of the cuprates. Neglecting the Ashkin- 
Teller interaction term present in this theory, the classi- 
cal model may be mapped onto a four-state clock model, 
with the basic variable being an angle param etrizing the 
four possible current loop oricntations.'SESIII] 

The quantum version of this model includes a kinetic 
energy term describing the quantum dynamics of the 
angle variables. Adding dissipation of angle differences 
as in the Caldeira-Leggett approach for Josephson junc- 
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tions, the model has been reported to exhibit local quan- 
tum criticality. Local quantum criticality in this context 
means that the model exhibits a fluctuation spectrum 
which only depends on frequency, but is independent of 
the wave vector P This essentially implies a dynamical 
critical exponent z — > cx3. A point which quite possibly 
is of importance in this context, is that while the starting 
point in Ref. 8 is a model with two Ising-like variables, 
the actual dissipative quantum model discussed is one 
with global U{1) symmetry. 

While the physical picture of fluctuating conflgurations 
of current loops suggests an identification of the angles 6 
and 9 + 2n, the presence of the clearly noncompact dissi- 
pation term makes this not entirely obvious. It is there- 
fore the intent of this work to investigate if the restric- 
tion of the variable domain influences the dynamical crit- 
ical exponent z, and thereby if it may have consequences 
for possible manifestations of local quantum criticality in 
similar models. Since it is still an open question exactly 
what the consequences are of how the variable domains 
are defined in dissipative quantum models, a numerical 
comparison of the compact and noncompact case is of 
general interest. We will therefore not restrict the inter- 
pretation of the model to Ising variables associated with 
loop currents, although the Z4 symmetry reflecting this 
starting point will be maintained. Moreover, due to the 
long-ranged interactions in the imaginary-time direction, 
the Monte Carlo computations are extremely demanding. 
Since we are interested in a proof of principle of the im- 
portance of compactness versus noncompactness, we will 
in this paper limit ourselves to a (1 + 1)D model. 

We will perform Monte Carlo simulations on two ver- 
sions of a dissipative Z4 model described in more detail in 
Sec. |lT] one with compact variables (i.e., a clock model) 
and one with noncompact variables. The simulation de- 
tails are described in Sec. after which we present the 
results, first for the noncompact case in Sec. |IV| then for 
the compact case in Sec. |V] 

Our main finding is that, although the critical scal- 
ing of space and imaginary time is equal for both cases, 
i.e., z = 1, there is a major difference in phase struc- 
ture. Whereas the compact model displays a conven- 
tional order-disorder phase transition, the noncompact 
model develops an intermediate phase characterized by 
power-law decay in spin correlations (quasi-long-range 
order) and a J7(l) symmetric distribution of the complex 
order parameter. The appearance of this intermediate 
phase is related to the fact that the kinetic energy term 
must be treated differently for the compact and noncom- 
pact cases, as we discuss in detail in the Appendix. 

It is well established that this kind of critical phase 
occurs in classical 2D Zg clock models and XY models 
with Zq anisotropyjE^HSr but only for larger values of q 
than we are considering. It is remarkable that the non- 
compact model presented in this paper exhibits a critical 
phase with emergent U{1) symmetry, when the dissipa- 
tionless starting point is a pure Z4 = Z2 x Z^ model 
(i.e., a double Ising model) with the angle variables re- 



stricted to four discrete values by a hard constraint. We 
will discuss this in more detail in Sec. VI after which we 
summarize our results in Sec. I VIII 



II. THE MODEL 

The starting point for our model is a chain of N^, quan- 
tum rotors, or equivalently planar spins, the alignment 
of which is described by a set of angle variables 
Although these variables could also be denoted as the 
phases of the quantum rotors, we will refer to them sim- 
ply as angles. Requiring that the spins satisfy sym- 
metry, the angles can be parametrized a& d — li^njAi with 
integer n, making our model similar to a four-state (or 
Z4) clock model. Being quantum spins, their dynamics is 
described by their evolution in imaginary time r, with Nr 
denoting the number of Trotter slices used to discretize 
the imaginary time dimension. The variables \dx^T\ a-re 
thus defined on the vertices of a (1 + 1)D quadratic lattice 
of size X Nr- 

In order to investigate if the restriction on the angle 
variable is relevant to the dynamical critical exponent 
z or not, we will consider two variants of this model, 
with the complete action for both stated below for later 
reference. In the compact (C) case, we restrict the 
parametrization variable n to just four values, so that 
the angle 9 is restricted to one primary interval, corre- 
sponding to the four primary states of the four-state clock 
model. In the noncompact (NC) case we have no such re- 
striction, and n can take any integer values. The general 
form of the action is 



S 



C,NC _ cC.NC 



— Sx -\- 5'dissj 



(1) 



where the kinetic energy for the compact and the non- 
compact case, respectively, is given by 



X=\ T=l 



X=\ T=l 

The spatial interaction defines a periodic potential 



(2) 



(3) 



5, = -ifg^cos(0,+i,,-0,,,), (4) 

X—l T—1 

and the dissipation term is defined according to 

2A.^2.^UJ sin^(^lr-r'l) ■ 

The bond variable or angle difference is written as 
Aflx.T = 9x+i,T — 9x.T- 
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Note that the only apparent difference between the 
compact and the noncompact model is the form of the 
kinetic energy term. When the angles are compact 
the short-range temporal interaction is given by a co- 
sine term, in contrast to noncompact angles for which a 
quadratic form of the kinetic term must be used. The 
reason for this difference can be traced to the fact that, 
whereas canonical conjugate variables of compact angles 
are discrete due to the 27r periodicity of the quantum 
wave functions, no such restriction applies when the an- 
gles are noncompact. From a qualitative point of view 
the two separate forms of the temporal interaction term 
is expected. Considering the imaginary time history of 
a single variable, it is clear that a cosine interaction in 
imaginary time will render the ground state of the non- 
compact model massively degenerate. A Trotter slice 
may be shifted by 2tt relative to the neighboring Trot- 
ter slices without any penalty in the action. However, a 
quadratic interaction term in the imaginary time direc- 
tion lifts this degeneracy and tends to localize the angle 
variables. 

There is nothing new about the derivation of these dif- 
ferent kinetic terms, but as the difference is crucial to 
the phase structure of our models and is also rarely dis- 
cussed in the literature, we include the derivation in the 
Appendix. In addition, in order to simulate the compact 
model we also need an appropriate reinterpretation of the 
dissipation term. We find it natural to postpone this to 
Sec. El 

The action is on a form identical to the model in Sec. 
Ill in Ref . [3] apart from the nature of the variables and 
the resultant treatment of the dissipation term. However, 
we still expect the scaling arguments presented in Ref. [31 
to be valid since no reference to the actual type of variable 
is used. The action in Fourier space may be written 

S ^{q' + cj^ + l^lq^ge-g, (6) 

neglecting any prefactors. Taking the limit q 0, uj 
we anticipate that the term ~ |w|q^ describing the 
dissipation is subdominant for all positive z. Accordingly, 
we expect at least naively that z ~ 1 for both compact 
and noncompact variables. This will be investigated in 
detail in our simulations, and we make no assumption of 
the veracity of naive scaling applied to this problem. 

III. DETAILS OF THE MONTE CARLO 
COMPUTATIONS 

When expanding the dissipative term, it becomes clear 
that it contributes both to ferromagnetic and antifer- 
romagnetic long range interactions. This renders the 
system intractable to the Luijten-Blotd^ extension of 
the Wolff cluster algorithnPSl which has been used with 
great success in systems with noncompeting interactions. 
Also, for the case of noncompact variables there does 
not exist a straightforward way of defining (pseudo)spin 
projections, a necessary point for the Wolff embedding 



technique.!^ Considerable progress has been made in con- 
structing new effective algorithms for long range inter- 
acting systems with extended variables .E^ti^ However, 
these algorithms are presently restricted to (0-1- 1)D sys- 
tems, and do not seem to generalize easily to Nx > 1.^^^ 
Furthermore, the basic degrees of freedom in these algo- 
rithms are the phase differences between two supercon- 
ducting grains in an array of Josephson junctions. Our 
aim is to investigate the ordering of the phases them- 
selves. Hence, the existing non-local algorithms may 
not be utilized. In the Monte Carlo simulations , we 
have therefore used a parallel tempering algorithnPfilll] 
in which several systems (typically 16 or 32) are simu- 
lated simultaneously at different coupling strengths. 

A Monte Carlo sweep corresponds to proposing a local 
update by the Metropolis-Hastings algorithm for every 
grid point in the system in a sequential way. For the 
case of noncompact variables the proposed new angles 
are generated by randomly choosing to increase or de- 
crease the value, then propagating the value by randomly 
choosing the increment on the interval { f , TT, ^,27r}. In 
the case of compact variables, a new angle value in the 
primary interval is randomly chosen. After a fixed num- 
ber of Monte Carlo sweeps (typically 3 — 10) a parallel 
tempering move is made. In this move, a swap of config- 
urations between two neighboring coupling values is pro- 
posed, and the swap is accepted with probability 'E.px 
given by 

"^^-\c-^ ifA>0. 

Here, A = k'{S[X; k']- S[X' ; k'])~- k{S[X; k]~ S[X' ; k]), 
where k is the coupling value varied, representing in our 
case K or a, and X represents the angle configuration. 
S indicates the term of the action proportional with the 
coupling parameter k. 

All Monte Carlo simulations were initiated with a ran- 
dom configuration. Depending on system sizes various 
numbers of sweeps were performed for each coupling 
value. For the phase transition separating the disordered 
state from the critical phase in the noncompact model 
5 — 10 X 10^ sweeps were made. Also, 1 — 5 x 10^ sweeps 
at each coupling value were discarded for equilibration. 
For the compact model and the second transition of the 
noncompact model as much as 30 x 10^ sweeps where 
performed and typically 5 x 10^ sweeps discarded. 

The Mersenne-Twisteil^^l random number generator 
was used in all simulations and the random number gen- 
erator on each CPU was independently seeded. It was 
confirmed that other random number generators yielded 
consistent results. We also make use of the Ferrenberg- 
Swendsen reweighting technique which enables us to 
continuously vary the coupling parameter after the sim- 
ulations have been performed. 



4 



IV. RESULTS: NONCOMPACT MODEL 

In this section wc consider the noncompact version of 
the dissipative Z4 model. Using Eqs. ([s]), Q and ([s]) we 
have the following action, 



a = 0.039 



a = 0.061 



S 



NC 



S: 



NC 



Sx + Sd 



(8) 



In contrast to the compact model, the angle variables 
are in this case not restricted to the primary interval. 
The variables are straightforwardly generalized to take 
the values 9 — 27rn/4, where n = 0, ±1, ±2, .... We seek 
to fix K and Kt and investigate how the system behaves 
under the influence of increasing dissipation strength con- 
trolled by the dimensionless parameter a. 

The kinetic coupling strength has been fixed to = 
0.4 for computational reasons, as this ensures that the 
simulations will be performed at convenient values of 
and Nt- We have performed simulations at four different 
spatial couphng constants K = 0.4, 0.5, 0.6 and 0.75. 
These choices are also made for computational conve- 
nience, as the limit of vanishing dissipation as well as 
the limit K ^ are both very computationally demand- 
ing. For all coupling values there is a disordered phase at 
low values of the dissipation strength. In this phase the 
noncompact angles exhibit wild fluctuations and conse- 
quently (e'^^'^) = 0. However, we also have (e'^^^-^) ^ 
in this phase, a trivial consequence of the cosine potential 
acting as an external field on the bond variables. The 
bond variables occasionally drift from one minimum of 
the extended cosine potential to another. As the dissipa- 
tion strength is increased, fluctuations in these variables 
are suppressed, and the system features two consecutive 
phase transitions separated by a critical phase. This in- 
termediate phase is characterized by power-law decay of 
spatiotemporal spin correlations on the form 



(9) 



The correlation functions for both spatial and imaginary 
time direction are shown in Fig. [T]for two different dis- 
sipation strengths both within the the critical phase. 

A very similar critical phase, as well as phase transi- 
tions associated with it, has recently enjoyed increased 
interest in various versions of classical clock models.'^SEH 
We will proceed under the assumption that a similar pic- 
ture is valid in our case. Indeed, simulations performed 
on a classical 2D six-state clock give qualitatively very 
similar results for all observables considered below, which 
supports the supposition that these two phenomena are 
related. 

Considering the complex order parameter of the sys- 
tem. 



(10) 



the intermediate critical phase can be identified by ob- 
serving the distribution of m in the complex plane. In 
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FIG. 1: (Color online) Correlation functions, Eq. ([9|, at 
two values of the dissipation strength q within the criti- 
cal phase for spatial coupling K = 0.75. System sizes are 
= 44,57,74, with optimal choices of A'',- at ai^'', see text. 
Top row: Correlation functions for the temporal direction. 
Bottom row: Correlation functions for the spatial direction. 



the disordered phase, the order parameter is a Gaus- 
sian peak centered at the origin. In the intermediate 
phase, quasi-long-range order develops in the complex 
order parameter, and so |77i| acquires a nonzero value as 
a finite-size effect. The order parameter is, however, free 
to rotate in the direction. This can be described as 
the vanishing of the excitation gap naively expected for 
discrete Zg models, or equivalently as an emergent U{1) 
symmetry.'-^^ This symmetry is broken at a larger value of 
the dissipation strength, when true long-range order is es- 
tablished when the magnetization selects one of the four 
well-defined directions in the complex plane originating 
with the underlying Z4 symmetry. Typical distributions 
of the complex order parameter in the three phases is 
shown in Fig. |2] 

Although not presented here, we have also confirmed 
that the susceptibility of the order parameter diverges 
over a finite interval of dissipation strengths, also a clear 
evidence of a critical phase. 

The phase transition between the disordered state and 
the intermediate critical phase at dissipation strength 



a = a. 
where 



(1) 



is detected by the Binder cumulant g — l — Q/3, 



Q 



{\m\') 

/|^|2\2 ■ 



(11) 



The brackets indicate ensemble averaging. The scaling at 
criticality of the Binder cumulant for anisotropic systems 
is given in terms of two independent scaling variableSj^^l 



Nx Nr 



(12) 



At a critical point the correlation length ^ diverges, and 
one should be able to observe data collapse of the Binder 



3R(m) 

(a)Two dimensional Gaussian distribution of the order 
parameter in the complex plane corresponding to the 
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(b)Intermediate critical phase exhibiting a 
finite-size-induced non- vanishing |m| that rotates in the 
<j> direction. The critical phase exists in a finite interval 
of dissipation strengths a^^^ < a = 0.04 < The 
remaining anisotropy is attributed to insufficient 
sampling.!^ 
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(c)The rotational symmetry of the intermediate critical 
phase is broken and long-range order is established as 
the order parameter relaxes into one of the four 
directions in the complex plane. The long-range ordered 
phase corresponds to the strong dissipation limit, 
a = 0.18 > oi^'. 

FIG. 2: (Color online) Evolution of the complex order param- 
eter when dissipation strength a is increased for K = 0.75 and 
system size A^^, = 74, Nt = 103 which corresponds to a near 
optimal aspect ratio at the phase transition at a « ac^\ The 
color scale indicates relative density of the distribution. 
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FIG. 3: The order parameter {\m\} and the anisotropy mea- 
sure of the noncompact Z4 model with K = 0.75 and 
system size Nx = 74, Nt = 103 ~ N*. This size represents a 
near optimal aspect ratio at q « a^^^ The two phase tran- 
sitions are indicated by arrows, note that the intermediate 



critical phase a 



(1) 



< a < a. 



(2) 



features a rotationally sym- 



metric order parameter distribution. 



cumulant as a function of N^/N^ for the correct value 
of z. The value of g{Nx,Nr) is independent of at 
the critical coupling, this may be used to align the plots 
of g as a function of N^. horizontally. The exponent z 
can then be found by optimal collapse of data onto a 
universal curve. The cumulant curves have a maximum 
at Nt — N* . At this temporal size, the system appears as 
isotropic as it can be, the anisotropic interactions taken 
into account. See Ref. |3]for a thorough discussion of this 
finite-size analysis. 

In the intermediate phase, the system is critical over a 
finite interval of dissipation strengths. According to the 
scaling Eq. ( 12 1 , curves of the Binder cumulant for in- 



creasing system sizes will therefore merge in this interval 
for Nx ooW^ For systems of finite sizes as considered 
here, the curves will however intersect close to the tran- 
sition instead, and we find ai^^ by inspecting the con- 
vergence of the crossing points. As discussed in Ref. |21 
the functional form of this convergence is unknown in our 
case (cf. also Sec. VI and Ref. [^^ . and all we can do 
is to report our best estimate for the — > 00 transition 
point. The uncertainty estimated accordingly is not in- 
significant, but the effective critical exponent z is found 
to not be very sensitive to this error in ac- 

By further increasing the dissipation strength, the ro- 
tational symmetry of the global order parameter is bro- 
(2) 

ken at a — ac . The Binder cumulant given by Eq. 
(11 1 will not pick up this transition because \m\ does not 



contain any information on the angular direction of the 
global magnetization. Ther efore , we consider an alterna- 
tive magnetization measur^^MZl 



m0 = (cos(40)). 



(13) 
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where 4> is the global phase as indicated by Eq. (10 1. 
This anisotropy measure vanishes when (j) is evenly dis- 
tributed and tends toward unity when the excitation gap 
opens and </> gets localized. We show in Fig. [3] both or- 
der parameters for the system = 74, Nr — 103 as a 
function of a. This N.,- corresponds to the nearest inte- 
ger N* Sit a K, a"c^. Actually, the optimal Nr decreases 
with increasing a, so the given system size does not rep- 
resent an optimally chosen aspect ratio for other dissipa- 
tion strengths. The rotational symmetry of the complex 
order parameter is clearly seen to be broken at a higher 
dissipation strength than the onset of the intermediate 
critical phase. 

Because (p measures a global rotation of the order pa- 
rameter, extremely long simulations is needed to explore 
the (f) space with a local update algorithm. This limits 
the efhciency of constructing a Binder cumulant from 

(2) 

and extracting from a universal point because this 
would involve calculating moments of a already statisti- 
cally compromised ensemble. To alleviate these difficul- 
ties, we instead make a scaling ansatz for the anisotropy 
measure itself. 
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FIG. 4: Phase diagram for the noncompact model, Eq. (|8| 
with Kt = 0.4. The dotted lines are guides to the eye. For 
fixed K the model features two consecutive phase transitions 
surrounding the intermediate critical phase (with quasi-long- 
range order). The simulation results (symbols along the dot- 
ted lines) are restricted to a region in coupling space amenable 
to simulations. 
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based on the fact that the naive scaling dimension of this 
magnetization measure is zero. Near criticality, we ex- 
pect to scale with system size in the same way as the 
Binder cumulant Eq. (12 1. Hence, we may calculate a 



dynamical critical exponent for this transition by exactly 
the same procedure as in Sec. |V]and Ref. 3j Again we 
expect a merging of curves as a — >■ ai^^ from above 
in the limit of large N^, but for the present system sizes 
we use the crossing points of curves to estimate ai^^ . 
In Fig. |4] we plot the resulting phase diagram in the 
a — K plane. The intermediate phase is evidently very 
wide also when compared to the uncertainty assigned to 
the transition line, and we feel confident that it is a gen- 
uine phase and not merely an effect of the admittedly 
moderate finite system sizes we are restricted to. 

We extract the dynamical critical exponent z along 
both of the critical lines ai^^ and ai^"* for all spatial 
coupling strengths. The data collapse of the Binder 
cumulant g at K = 0.75 and a = 0.030 w ai^^ is 
shown in Fig. |5] Increasing the dissipation strength fur- 
ther brings the system to the second phase transition at 

(2) 

a — 0.125 Q!c , the collapse of at this point is 
shown in Fig. |6] 

In Table |T] we present the numerical estimates of the 
dynamical critical exponent. The values of z are obtained 
using the scaling relation N* — aN^, with uncertainties 
based on a bootstrap analysis. These uncertainties also 
include the uncertainty in ac- Within the accuracy of the 
simulations, the value of the critical exponent is z = 1 
for all the coupling values at both phase transitions (al- 
though precise results are harder to obtain for the sec- 
ond). This is in accordance with the scaling argument 
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FIG. 5: (Color online) Data collapse of the Binder cumulant, 
(7=1 — Q/3, with Q given by Eq. (Ill, for the noncompact 
^4 model &t K = 0.75 and a = 0.030 ^ a*^' with z*^' = 1. 
Inset: Intersection of the Binder cumulant as a function of 
dissipation strength. 



presented in Sec. ITTl 



RESULTS: COMPACT MODEL 



We now turn to the compact version of the dissipative 



Za model, 



(15) 



where the three terms are given by Eqs. (|2|, (Q and (|5|, 
respectively. Note that we now use a kinetic term 5*^ hav- 
ing the same cosine-form as the spatial interaction term 
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FIG. 6: (Color online) Data collapse of the anisotropy mea- 
sure m^, Eq. ( 13 1, for the noncompact Z4 model at 7f = 0.75 
and Q = 0.125 ~ OL-if^ with z'^' = 1. The actual uncer- 
tainties are probably larger than indicated by the error bars 
for reasons discussed in the text. Inset: Intersection of the 
anisotropy measure as a function of dissipation strength. 



TABLE I: Numerical estimates for critical coupling and crit- 
ical exponents z^^-*''^' for the two phase transitions ac^'''^' of 
the noncompact model. 



K 




^(1) 


(2) 
etc 




0.75 


0.030(2) 


0.99(1) 


0.125(2) 


1.01(2) 


0.6 


0.042(2) 


0.99(2) 


0.190(3) 


0.96(3) 


0.5 


0.053(2) 


1.02(2) 


0.238(3) 


0.97(3) 


0.4 


0.068(4) 


0.97(3) 


0.287(5) 


0.99(4) 



Sx ■ Regarding the use of the same dissipation term S'diss 
as in the noncompact case, one may argue that adding 
a Caldeira-Leggett term for the angle differences AS is a 
rather artificial way to model dissipation for a compact 
clock model in the first place, since its variance under 
27r translations of Q implicitly assumes noncompact vari- 
ables. However, adding exactly such a dissipation term is 
crucial for the demonstration of local quantum criticality 
in a similar modePthat is not obviously noncompact. 
Therefore, our motivation for the comparative study in 
the present section of a compactified version of the ac- 
tion ([8| is to investigate whether an equivalent dissipa- 
tion term for compact variables gives the model the same 
critical properties as reported for noncompact variables 
in the previous section, and thus whether the compact- 
ness of the variables as such is essential. Constructing an 
appropriate compactified version of the dissipative model 
does, however, require a reinterpretation of the variables 
in the Caldeira-Leggett term, so we will begin with a 
careful discussion of how we should treat this term in 
our simulations. 



We first impose the following restriction on the in- 
terpretation of the compactified dissipation term: The 
term as a whole should be invariant under translations 
B — >■ 6* -I- 27r, since these two states are indistinguish- 
able. As a corollary, any configurations that are phys- 
ically indistinguishable when the angles are restricted 
to four values Q G {— tt, — 7r/2, 0, 7r/2} (or any equiva- 
lent parametrization) should give the same contribution 
to the dissipation term. Consequently, we cannot sim- 
ply simulate the model with the dissipation term (|5| as 
it stands, because the angle differences AO^^t now only 
make physical sense modulo 27r. We therefore have to 
bring AOx.t back to the primary interval [— 7r,7r), as is 
well known for phase differences in superconducting sys- 
tems without dissipation and other realizations of the 
(compact) XY model. Furthermore, we also choose to 
do the same for the difference between the two (com- 
pactified) A9x^T terms in Eq. ([s]), as the alternative 
would result in different Boltzmann factors being asso- 
ciated with physically equivalent situations. Our proce- 
dure then is equivalent to requiring that the entire differ- 
ence A0x^r — AOx^t' should be restricted to the primary 
interval [— 7r,7r), i.e., treating the dissipation term as a 
27r-periodic function. 

The details of the Monte Carlo simulations are de- 
scribed in |III| also for the compact model. The only dif- 
ference that may be of any consequence is that we found 
it more convenient to vary the spatial coupling while fix- 
ing the dissipation strength in this case, but we have 
checked that the direction in coupling space taken by the 
simulations has no impact on the result. 

The dissipationless (a — 0) four-state clock model 
is completely isomorphic to the Ising model with in- 
teraction K/2. Thus, we may employ the criterion 
sinh(Xc) sinh(ifT-) = 1 in order to calculate Kc for a fixed 
value of K^. The temporal coupling parameter is fixed 
at Kr = — In (tanh i) « 0.7719 such that Kc — \ when 
the dissipation is tuned to zero. 

The most striking difference we found when compacti- 
fying the angles is that the intermediate phase with quasi- 
long-range order vanishes. This means that one has only 
a single disorder-order phase transition, as is the result 
one would usually expect for any model with sym- 
metry. We have verified that the symmetry and the 
apparent U{\) symmetry of the complex order parame- 
ter (in the disordered phase) are spontaneously broken 
simultaneously at a single critical point. This is found 
by observing that the inflection points of magnetization 
curves for m and coincide asymptotically, in contrast 
to the curves shown in Fig. [3] for the noncompact case. 

The phase diagram for the compact Z4 model with 
bond dissipation is shown in Fig. [7] It differs consider- 
ably from that of its noncompact counterpart, not only 
in the evident absence of any intermediate critical phase, 
but also in that the limit a — >■ is well-behaved. Here, 
the model is reduced to two uncoupled 2D Ising models, 
for which exact results are known and simulations are 
straightforward. In the limit of — >■ the simulations 
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FIG. 7: Phase diagram for the compact model of Eq. (151 
with Kr = — ln(tanh |), the dotted line indicating a critical 
line separating the disordered phase from a phase with long- 
range order. The line is not drawn beyond a = 0.4 because 
of increasing uncertainties. 



TABLE II: Critical coupling Kc and dynamical critical expo- 
nent z for different values of the dissipation strength a for the 
compact model. 



a 




z 


0.05 


0.8303(4) 


1.02(2) 


0.1 


0.6753(7) 


0.99(2) 


0.2 


0.414(2) 


0.99(2) 



are on the other hand very difficult for the same reasons 
as those investigated by us in a similar model in Ref. 
[21 Therefore, we have not strived to extend the phase 
diagram all the way down to the a. axis in this work. 
Due to the qualitative difference in the kinetic terms for 
the compact and noncompact model, it is not possible 
to make quantitative comparison between the position of 
the phase transition line in Fig. [7] and the two phase 
transition lines in Fig. [4] 

Turning next to the nature of the critical line in the 
phase diagram, we show in Fig. [8] and Table |TT] results 
for the three points along the line for which we made the 
most effort to extract the dynamical critical exponent. 
These points are chosen so that the relative influence of 
the dissipation term should be qualitatively comparable 
with that for the points (ai^\ K) chosen for the first tran- 
sition of the noncompact model. As for the noncompact 
model here and the Ising model with bond dissipation 
studied in Ref. [31 there is no significant variation in 
the dynamical critical exponent from the expected value 
z although the tendency to greater finite-size effects 
for increasing a remains for both the compact and the 
noncompact model. 



FIG. 8: Finite-size analysis of the maximum A'^^ of the cu- 
mulant curves as a function of spatial system size Aj, used 
to obtain the dynamical critical exponent z for the compact 
model. The dashed lines show the power- law fits, cf. Table [ill 
for the results. 



VI. DISCUSSION 

The models discussed in this paper are in some sense 
generalizations of the Ising spin system with bond dissi- 
pation discussed in Ref. T For the compact Z4 model the 
modifications come from the increase of the number of 
states from g = 2 to 5 = 4, while one for the noncompact 
model adds an additional extension of the configuration 
space. The phase diagram for the compact model is very 
much like that observed for the dissipative Ising modeljl^ 
both featuring a single order-disorder phase transition 
line. The noncompact model on the other hand exhibits 
much richer physics in the sense that it presents, for fixed 
K and , two phase transitions surrounding an interme- 
diate critical phase with power-law decaying spin corre- 
lations and emergent f7(l) symmetry. The most pressing 
question then pertains to the occurrence of this phase: 
Why is the discrete structure of the angle variables ren- 
dered irrelevant in a region of parameter space for our Z4 
model, when such behavior is previously known to occur 
only in Zq models with (7 > 4? Even our compactified 
model differs from a pure Z4 = Z2 x Z2 clock model, 
since the dissipation term couples the two underlying Z2 
models in a nontrivial way. Such models can no longer be 
expected a 'priori to behave as an Ising model, and there 
is in principle no reason why they may not even present 
intermediate phases. The absence of such a phase in our 
compact model does however indicate that we must turn 
to the other obvious difference between our model and a 
Z4 clock model, namely that the variables in our noncom- 
pact model are free to drift outside the primary interval. 
Somehow, this added degree of freedom is enough to close 
the excitation gap. 

As observed in Fig. 2(b)l the underlying Z4 symme- 
try stemming from the discreteness of the variables is 
irrelevant in the intermediate phase. Consequently, the 
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system displays an effective continuous symmetry. Since 
z = 1, the effective long- wavelength low-energy propa- 
gator is on a Gaussian form l/(w^ -|- q^). In addition, 
the system is effectively two dimensional due to z = 1. 
In two dimensions, Gaussian fluctuations are sufficient 
to induce a critical phase given a continuous symmetry. 
This is analogous to the mechanism producing a criti- 
cal phase in the classical 2D XY model with an Zg>4 
anisotrop}^ [soft constraint with underlying U(Vj sym- 
metry], and also for classical ^q>4 clock models^ (hard 
constraint). The difference in our case is that the under- 
lying symmetry is Zq^^. 

To comment further on the origin of the critical phase, 
it appears that the quadratic form of the kinetic energy in 
the problem is essential for observing it. This quadratic 
short-range interaction term in imaginary time facilitates 
Gaussian fluctuations. Were we to use a cosine-like form 
of this term for noncompact variables (as one does for 
compact variables), this intermediate phase would not 
be found. The kinetic energy term is bounded from be- 
low, but not from above. Upon entering the intermediate 
phase from the ordered side, this term tends to suppress 
strong 9 fluctuations, much more so than a kinetic term 
which is bounded from below anrf above, such as a cosine- 
like term. Only at even lower values of the dissipation are 
the excitation energies of larger 9 fluctuations so low that 
wild 9 fluctuations are possible due to the boundedness 
of the spatial coupling. At this point, the system dis- 
orders completely. If the quadratic kinetic energy term 
is replaced by a cosine-like term, wild 9 fluctuations are 
facilitated precisely at the critical point where the Z^- 
symmetry becomes irrelevant, and the system disorders 
directly from the Z4-ordered state. Hence, for a compact 
model there will only be one phase transition separat- 
ing the Z4-ordered state from the completely disordered 
phase. 

We now comment on the critical scaling between space 
and imaginary time in the models we have studied. For 
the compact case one has the conventional case of a crit- 
ical line along which the correlation length diverges as 
( ^ \K - Kc]-" in space and \K - Kd^^" in 

imaginary time, with z appearing to remain equal to 
unity along the line. This picture is no longer valid 
in the noncompact ^ and are formally infl- 

nite in the entire intermediate critical phase, and z can 
not be defined from the anisotropy of their divergence 
in this region. Furthermore, supposing that the inter- 
mediate phase shares qualities with the corresponding 
phase in classical Zq^^ models, the correlation lengths 
can be expected to diverge exponentially as this critical 
phase is approached from either side as for the Kosterlitz- 
Thouless (KT) transition, and not as a power law as for 
conventional critical points. However, as long as the cor- 
relation length does diverge, and this divergence is expo- 
nential both in space and imaginary time, the dynamical 
critical exponent is still well defined through ~ 

Therefore, our finite size analysis is valid as a — > ai^"^ 



and a — 7> ac irrespective of whether these points turns 
out to possess KT criticality or not. At both phase tran- 
sitions we have z — I, signaling equally strong divergence 
of correlation lengths in space and imaginary time. 

To infer from simulations on finite systems that the 
correlation length in fact diverges exponentially is ex- 
ceedingly difficult jf^SlEil and we have not attempted to 
determine the exact nature of the phase transitions, but 
leave this an open question. The phase transitions (one or 
both) may be in the KT universality class, or it may be- 
long to a class of related topological phase transitions .f^Sl 
This identification of the exact universa lity clas s is con- 
troversial even for classical clock models .^^^^^Ml 

If we generalize the noncompact action in Sec. |IV| by 
redefining the phase space such that the variable can take 
on all real values, Eq. (|8| may represent the action for a 
one-dimensional array of Josephson junctions.'^ Recent 
theoretical worlP^ESl report that such systems may dis- 
play local quantum criticality, in the sense that the spa- 
tial coupling renormalizes to zero at the quantum phase 
transition so that the behavior is essentially (0 -f 1) di- 
mensional. This suggests that local quantum criticality 
need not be restricted to (2 + 1)D models such as the 
one presented in Ref. |S1 but that similar unconventional 
criticality may be found in (1 -|- 1)D as well. Although 
it should be remembered that our (1 + 1)D model has 
discrete angle variables, our simulations do not show any 
traces of local critical behavior, in the sense that the scal- 
ing of Binder cumulants do not give z 3> 1 . 

Strictly speaking, the dynamical critical exponent is 
not well defined inside the intermediate phase, and the 
isotropic behavior is instead maintained by the decay ex- 
ponents for the power-law spin correlation functions in 
space and time being equal. Nevertheless, for finite 
one may still assume the scaling relation N* = aN^ and 
use the ordinary procedure to extract the (effective) ex- 
ponent z as long as the system is critical, which yields 
z « 1 in the entire intermediate phase. We may then 
inspect how the nonuniversal prefactor a changes as a re- 
flection of the anisotropy of the interaction in time and 
space. In the noncompact model it is possible to investi- 
gate the development of a at constant K^/K and varying 
a without leaving the critical region. We find that a de- 
creases for increasing a, indicating that the dissipation 
term contributes to making the temporal dimension less 
ordered than the spatial one. This is also in contrast 
with a tendency toward (0 -I- 1)D behavior when increas- 
ing the dissipation strength, as suggested in the models 
mentioned above. 



VII. CONCLUSIONS 

We have performed Monte Carlo simulations on two 
distinct Z4-symmetric dissipative lattice models. In one 
model the phase variables are only defined on the inter- 
val [0,27r), while the other model has no restrictions on 
the variables. The different domains of the variables have 
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implications for the short range interaction term in imagi- 
nary time, which again leads to essential differences in the 
behavior of the two models. The compact model features 
only one phase transition in which the Z4 symmetry is 
spontaneously broken. On the other hand, the noncom- 
pact model displays three phases, namely a disordered 
phase with exponentially decaying spin correlations, an 
intermediate critical phase with quasi-long-range order, 
and finally a long-range ordered phase. 

Along the phase-transition line of the compact model, 
we find the dynamic critical exponent z = 1, independent 
of the dissipation strength. In the noncompact model, we 
find the value z = 1 for both phase transitions and the 
power-law decay exponents for space and imaginary time 
are equal in the entire phase exhibiting quasi-long-range 
order. 

We have shown that the issue of compactness ver- 
sus noncompactness of the fundamental variables of the 
Z4 models have important ramifications for their long- 
distance, low-energy physics. 
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Appendix A: Quantum-to-classical mapping for 
compact and noncompact variables 

In this appendix we will outline the quantum-to- 
classical mapping for a quantum rotor model and show 
how the kinetic term in the resulting classical model de- 
pends on whether the variables are interpreted as com- 
pact or noncompact. We will first reproduce the deriva- 
tion as given in Refs. I36I37I for the case of compact 
variables, after which we will generalize and reinterpret 
it for the noncompact case. Although there is nothing 
novel about this derivation, the form of the kinetic term 
often seems to be taken for granted in the literature, and 
a correct interpretation of the classical action in the non- 
compact case is crucial for our results. As a starting point 
we take the (dissipationless) Hamiltonian Hq ^ T+U for 
a spatially extended system of particles, each moving on 
a ring. The kinetic energy of the rotors is given by 



rotors, 



r = - 



1 



92 



21 ^ del 



(Al) 



where / is some inertia parameter. The (periodic) po- 
tential energy is given by Josephson-like coupling of the 



U = -i\:^cos( 



(A2) 



with K being the coupling strength. Here we have used 
the angle representation where we for simplicity let be a 
continuous variable, and B is the corresponding operator. 
Characteristic of a rotor model is the invariance of the 
system upon translations of the angle Q — >■ -I- 27r. The 
eigenfunctions describing the system should therefore be 
27r-periodic, a requirement which immediately yields dis- 
cretized angular momenta and energy levels. 

The partition function of the rotor system may be 
given by 



Z = Tr (e-^(^+^)) 



(A3) 



We let fcs — 1 such that /3 equals inverse temperature. 
The trace may be evaluated by introducing a path inte- 
gral over M time slices between r = and t = f3, with 
the width of the time slices given by At = P/M . For 
every time step indexed by t, we insert a complete set of 
states, 

lim veW {e{T + l)\e-^^^e-^^^\e{T)). 

•' r=0 

(A4) 

Here, |6'(r)) is an angular eigenstate of all rotors with 
Trotter index t. Since \0{t)) is an eigenstate of 9 we get 



'\e{r)) = \e{T))e' 



(A5) 



A general matrix element describing the kinetic energy 
is given by 



T.,r = {e.{r+l)\e-'^^'%{T)) 



(A6) 



Next, for each r we insert a complete set of eigenstates 
of the kinetic energy \nx(T)). Because 6 and n are con- 
jugate variables, we have the identity {nx{T)\Ox{T)) — 
exp [— irij;_,-0a;_7-]. Inserting this, we get the general form 
of the matrix clement for the kinetic energy 



■^e 2J 



(A7) 



Using the Poisson summation formula, we may write the 
summation over integer valued angular momenta in Eq. 
( A7 1 as an integral over the continuous field h at the cost 



of introducing another summation variable m: 

00 „ 

T^^r^ E / dne'"(^--+i-^=--'-57'^™'e^"'"" (A8) 



m— — 00 
00 



?n — — 00 
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where K^^ 



I 

At ' 



and C 



is a constant prefactor over the eigenstates in Eq. (|A7]) has to be replaced by an 



which is henceforth dropped from the expressions. The 
last approximation of Eq. ( A8 1 is the Villain approxima- 
tion of the cosine function, which is known not to alter 
the universality class of the phase transition. 

Reintroducing the matrix elements to the partition 
function and renaming Arif — > K , we get 



Z = j VOe^^ cos(e,,^+i-6 



(A9) 



i.e., an anisotropic XY model in (1 + 1) dimensions. Note 
however, that we were able to cast the kinetic energy ma- 
trix element into the form of a sequence of Gaussians 
because the angular momentum eigenvalues where re- 
stricted to integer values. This is only the case when the 
canonical conjugate variable is restricted to a [0, 27r) 
interval. In other words, the partition function given in 
Eq. (A9l reflects the interpretation of Eqs. (lAll) and 



(A2) in terms of rotors. 



Equations (Al) and (A2| may also describe particles 



moving in an extended potential, in which case the state 
of the system after a 27r translation is distinguishable 
from the state prior to the translation. Introducing dis- 
sipation to this system by coupling AO to a bosonic bath 
explicitly breaks the periodicity of the quantum Hamilto- 
nian, and consequently the variable 9 should be treated as 
an extended variable from the outset. This necessitates 
a modification of the above procedure as the summation 



integral over a continuum of momentum states. Then, 
the kinetic energy matrix element instead becomes 



(AlO) 



where a constant factor has been ignored. Inserting this 
expression into the kinetic part of the partition function 
yields 



Zr = lim / V0e 



2 L^T = 



(All) 



This continuum expression for the action is the one con- 
ventionally stated in the literature both for compact 
and noncompact variables. However, it is always im- 
plicit that the imaginary time dimension is discrete by 
construction,'^ and for most numerical computations it 
has to be treated as such in any case. One then has to 
choose one of two alternative discretizations of the short- 
range interaction in the imaginary time direction, de- 
pending on the interpretation of the system and the com- 
pactness of the variables. As shown above, the cosine- like 
term of Eq. (|2| is the natural discretization for compact 
variables, whereas the quadratic term used in Eq. (|3]) is 
associated naturally to noncompact variables. 
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